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This thesis investigates the use of the step frequency waveform, its design and analysis 
using the ambiguity function. The step frequency waveform consists of a series of N pulses 
each with a pulse width of x, and whose frequency is increased from pulse to pulse in steps 
of Af. A design procedure for detection of small targets with a surface (land or sea) based 
step frequency radar employing a high pulse repetition frequency (PRF) waveform is 
developed. The proposed method determines the waveform parameters for given radar 
specifications. A simple graphical implementation as well as a computer implementation are 
presented. The theoretical dimensions of the step frequency waveform are defined and verified 
for some waveforms with parameters similar to the waveform of interest. Finally, the 
ambiguity function is used to analyze the step frequency waveform. 
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I. INTRODUCTION 



A. BACKGROUND 

Over the past few years the applications for imaging or multidimensional high 
resolution radar systems have expanded greatly. This has been primarily due to the 
improvements in both the waveform generation and the signal processing needed to support 
the radar imaging technology. Higher resolution measurements give more information about 
the target, but also require more detailed measurements of the target, thus translating into 
increased radar system requirements and complexity. The theory of high resolution radar has 
been around for some time, but extensive applications of the technology have appeared only 
recently due to new developments in technology. Recent applications of high resolution 
techniques have appeared in a variety of areas, including: 

• Modeling and analysis or radar targets; 

• Diagnostic methods for target radar cross section reduction; 

• Development of target detection, discrimination, recognition, and 
classification techniques; 

• Improved radar detectability performance; 

• Diagnostic methods to support the medical and physical sciences. 

High range resolution (HRR) in radar has many good uses. Among these a few are 
listed below: 

• Enhancement of target-to-clutter ratio; 

• Accurate range measurement; 

• Resolution of multiple targets - raid count; 
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• Target range profile and size; 

• Target classification; 

• Recognition and avoidance of deception ECM; 

• Measurement of range rate without doppler. 

HRR is achieved by using a large bandwidth. In ultrawideband (UWB) waveforms, large 
bandwidth is obtained by decreasing the pulse width. However, when the pulse width is too 
narrow, sampling becomes problematic and is limited by existing analog-to-digital (A/D) 
technology. In conventional pulse compression (PC), a large bandwidth is obtained by 
modulation inside the pulse. In both approaches, the instantaneous bandwidth is large which 
requires a high A/D sampling rate. In the step frequency waveform, the bandwidth is 
increased sequentially (instead of instantaneously as in the previous cases) over many pulses 
within a burst by stepping the frequency of each pulse. By keeping the instantaneous 
bandwidth low, the requirements for A/D sampling rates are lowered considerably. The A/D 
sampling rate is a bottleneck in HRR radar design. Thus , the step frequency seems to be a 
good waveform for HRR. 

A step frequency radar operates by transmitting a waveform composed of a number 
of increasing stepped frequencies. A coherent I-Q detector measures both the in-phase and 
in-quadrature components of the received signal at each of the frequency steps. The phase 
information included in this complex sequence allows for accurate distance measurements. 
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B. OVERVIEW 



This thesis discusses the step frequency waveform design and its analysis using the 
ambiguity function. In Chapter II the fundamental principles of the step frequency radar are 
explained. The step frequency waveform is defined and the waveform processing is explained. 
Also in this chapter the implementation of the radar system is described. In Chapter HI the 
waveform design method is developed. Both graphical and computer implementations are 
shown and used to compute the waveform parameters. The ambiguity function is introduced 
in Chapter IV. Some important examples of the ambiguity function are shown. The theoretical 
dimensions of the ambiguity diagram are defined and then verified with a few cases of 
interest. These theoretical dimensions are then used to define the ambiguity diagram of a 
waveform with a large number of pulses, which for this reason is difficult to compute and plot 
with the desired resolution. Finally, in Chapter V the conclusions of this thesis are presented. 
The MATLAB codes used are given in Appendix B. 
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II. STEP FREQUENCY RADAR - OPERATING PRINCIPLE 



A. STEP FREQUENCY WAVEFORM 

The step frequency wavefonm consists of a series of N pulses each with a pulse width 
of x, and whose frequency is increased from pulse to pulse in steps of Af, as shown in Figure 
1 . The frequency of the w* transmitted pulse is given by 

fn "/o + («-l)¥, (2.1) 



where 



/ 0 = nominal carrier frequency; 
A /= frequency step size; 
«=1,2,3,...,N. 



f 0 - Nominal carrier frequency 
A f~ Frequency step size 




PRI 



Figure 1 . Step Frequency Waveform. 
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The total effective bandwidth of the step frequency waveform is given by 



B eff = W- 



( 2 . 2 ) 



The range resolution for this waveform is not the conventional value of ct/2. The 
synthetic range resolution obtained by coherently processing a series of N returns from the 
target is given by 



where AR is the processed range resolution in meters and c is the speed of light (3xl0 8 m/sec). 
As can be seen from Equation 2.3, the range resolution can be improved by either increasing 
the number of pulses or the frequency step size. 

The coherent processing interval (CPI) is the total time duration in which returns 
from a target are collected for signal processing. Thus the CPI is just the product of the 
number of pulses N, and the pulse repetition interval (PRI); 



The coherent processing interval is also known as the burst time. 

B. SYSTEM DESCRIPTION AND IMPLEMENTATION 

The implementation of the step frequency radar is similar to the coherent Pulse- 
Doppler radar. The system implementation can be observed in Figure 2. The core of the 
system is a coherent step frequency synthesizer with an output frequency called which 




(2.3) 



CPI = N(PRI). 



(2.4) 
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/« ~ f STAIO + fsrs + fcoao 



STALO 



Up-con vertcr 
Mixer 



Amplifier Circulator 




Figure 2. Block Diagram of a Step Frequency Radar. 
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is stepped from pulse to pulse by a fixed step size, A /. Therefore, we can write 

fsyn = (2.5) 

The frequency of a local coherent oscillator (COHO ),f coho , is first mixed with the 
output frequency of the step frequency synthesizer,/^ . Then the sum frequency is converted 
to the final transmitted frequency by mixing it with the frequency of the stable local oscillator 
(STALO), f stalo . The sum frequency of this second mixer is then pulse modulated and 
amplified before the transmission occurs. Each transmitted pulse will have a carrier frequency 
comprised of three different components: the fixed IF frequency, f coho> the fixed RF frequency, 
f 5talo , and the variable frequency of the stepped frequency synthesizer,/^ . Thus, the carrier 
frequency of the pulse is given by the following equation: 

fn = f stalo + f :oho + f syn = fo + ("“W- (2.6) 

In the receiver side, the return signal is first down converted by mixing it with f stalo . 
The resulting signal is further down converted to intermediate frequency (IF) by mixing with 
. The output signal obtained after the second mixer is the IF signal. The signal is then 
divided into two different channels, the in-phase channel (I), and the quadrature channel (Q). 
Both, the in-phase and the quadrature signals are mixed with the output signal of the 
coherent oscillator. The first is mixed directly and the second is mixed with the COHO signal 
after this last one has been shifted 90 degrees. The I and Q components, both in the video 
frequency range, are sampled by an A/D converter at a rate equal to the inverse of the pulse 
width. 
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Each complex sample is usually called range bin , thus range bin may be defined as a 
memory location for temporary storage of successive pairs of numbers representing the I and 
Q samples of the radar return received at a given point in the interpulse period. A separate 
bin therefore must be provided for each sampling interval (range gate). To the extent that 
range is unambiguous, the numbers stored in any one bin represent successive returns from 
a single range increment, hence the name range bin. Because of the correspondence of the 
range bins to the sampling intervals (when A/D conversion follows I and Q detection), range 
bin has come to be used synonymously with sampling interval as well as range gate. 

Sampled returns from incoming pulses are stored im memory until all the pulses within 
the same burst have been received and can be processed. Figure 3 describes the organization 
of sample storage in memory. Complex samples for each range bin are transformed by taking 
the FFT of samples in each range bin, resulting in a high resolution range profile (HRR 
profile). 

One important advantage of the step frequency radar compared with other radars 
which also use wideband waveforms, is that in the former the narrow instantaneous bandwidth 
does not require a high analog to digital sampling rate which can be a major limitation in the 
system design. On the other hand, the step frequency radar may require more complex signal 
processing. 
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n 



f„+Af 



f 0 +2Af 



f„+3Af 



f n +4Af 



f n +5Af 



N samples 



1 

1 



i 



1 



1 



t 



N 



f 0 +(N-l)Af 



r~L 



Range cells 






Range bin 



Figure 3. N samples from a target are collected and stored in memory for 
processing. 
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C. WAVEFORM PROCESSING 



The radar return signal at the output of the coherent detector can be written as 

s(n) = I(n)+jQ(n) = A n exp(j<$> n ) , (2.7) 

where 

A„ = the magnitude of the n th pulse; 

4>„ = the phase of the pulse. 

A„ represents an amplitude factor which depends on the transmitted power, the size of the 
target, its dielectric constant (relative to that of the medium where it is embedded), and the 
propagation losses in the medium. The signal components are shown in Figure 4. 





A, 


a > 


k I (n) = A„ cos 4> n 

X i 




\ Q( n ) = sin <j>„ 

/X \ ‘Pn w , 






In 



Figure 4. Vector Representation of the Signal Components. 
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The phase of the «* pulse can be written as 



<f>„ = InfJr = 2icf„ — (2.8) 

c 

where /„ is the carrier frequency of the w* pulse and t r is the round trip time from the target 
measured by the pulse. This time is equal to 2 R/c , with R n being the range between the 
target and the radar. For a target with constant velocity toward the radar, the range R„ is 
given by 

R „ = R o-( n ~ l ) vT ’ (2.9) 



where 

R 0 = initial target range; 

T = pulse repetition interval; 
v = relative radial velocity between target and radar. 

The overall expression for the baseband return signal can now be determined by 
substituting the instantaneous frequency from Equation 2.1, and the target range from 
Equation 2.9 in Equation 2.8, and then substituting the resulting expression in Equation 2.7. 
The resulting equation is as follows: 



s(ii) = A n exp 
= A „ ex P 



J~f n R n 
c 



j— (/„♦(»- 1)A/) (V(n-l)v7) 



( 2 . 10 ) 
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The set of N sampled baseband signals s(n) from a given range bin is converted to the 



time domain range profile by taking the FFT; 



H, = — Y \s(n) expf w/7/? 0<k<N-\. 
k N T N ) 



( 2 . 11 ) 
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III. STEP FREQUENCY WAVEFORM DESIGN 



A radar can be classified as low PRF, high PRF or medium PRF. A low PRF is, by 
definition, one for which the maximum range the radar is designed to handle lies in the first 
range zone - the zone from which first-time-around echoes are received. Range for low PRF 
is unambiguous and doppler frequency is ambiguous. A high PRF is one for which the 
observed doppler frequencies of all significant targets are unambiguous and range is 
ambiguous. Finally, a medium PRF is one for which neither of these conditions is satisfied. 
Both range and doppler frequency are ambiguous. In this thesis the interest is the high PRF 
radar. 

The purpose of this chapter is to develop a design procedure for detection of small 

targets with a surface (land or sea) based step frequency radar which employs a high PRF 

waveform. The disadvantage of high PRF is that the range ambiguity has to be resolved. 

However, several advantages result from high PRF. The unambiguous doppler leads to 

advantages such as elimination of clutter and absence or minimization of the blind doppler 

problem. It also improves the signal-to-noise (S/N) ratio by integration of larger number of 

pulses. For a constant frequency high PRF waveform, multiple time around clutter (from 

ambiguous regions due to different pulses) adds up and thus decreases the signal-to-clutter 

(S/C) ratio and increases the dynamic range of the receiver and associated A/D. Though the 

S/C ratio can be improved by doppler filtering, receiver and A/D will still require a large 

dynamic range as a large magnitude of folded clutter has to be accommodated in the receiver 

and A/D. The step frequency waveform reduces the folded clutter and thus the dynamic range 

by decreasing the number of range ambiguities. Range ambiguities in the step frequency 
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waveform depend upon the number of pulses which will pass through the EF filter 
simultaneously. This number is given by l/(xA f). By proper choice of the product xA /, range 
ambiguities and thus clutter and dynamic range can be reduced. However, this advantage 
comes about at the cost of limiting the range coverage. This feature of the step frequency 
waveform provides a good tradeoff as high resolution radars are often limited by A/D 
technology. The waveform design involves the selection of PRF (/ r ), number of pulses (N), 
frequency step (A/) and pulse width (x). It is assumed that the following radar performance 
related parameters are available: range resolution (AR), nominal carrier frequency (f 0 ), time- 
on-target (tot), minimum and maximum radial velocity (v^ and ) and the maximum 
number of pulses (N^J. The design of constant frequency waveforms is relatively simple and 
well documented. However, the design of step frequency waveforms is complex due to the 
conflicting effect of waveform parameters on various desired performance goals. Thus, in this 
chapter effort is made to systematize the design process for the high PRF step frequency 
waveform. The main point behind the design of high PRF waveforms is to keep the target in 
clutter free area and the target doppler unambiguous. This involves selection of PRF based 
on minimum and maximum target velocities and other parameters. PRF should be large 
enough that target at maximum expected velocity stays unambiguous in doppler. The other 
constraint is based on the fact that target doppler from minimum target velocity exceeds the 
clutter doppler. These constrains along with resolution requirements lead to the high PRF 
waveform design. 

In this chapter, the method used to compute the desired parameters for the high PRF 
design case will be explained, but first two parameters will be defined and the expression for 
constraints on PRF will be derived. 
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1 . 



Definition of two constants 



A dimensionless parameter P that relates the radial target motion over the time 
duration of the step frequency waveform to the resolution of the profile can be defined as: 



p - 

A R ’ 



(3.1) 



where: 

v = the radial velocity between the radar and the target (closing velocity); 
NT = the duration of the step frequency waveform; 

AR = the processed range resolution. 



This is the ratio of how far the target moved with respect to the range resolution cell during 
the integration time. For all velocities with magnitude greater than zero, this represents a 
mismatch in the DFT (which acts like a matched filter), resulting in attenuation and dispersion 
of the HRR profile. In addition to the attenuation and dispersion that will result in a reduction 
in range resolution and S/N ratio, the uncompensated radial velocity will shift the HRR 
profiles by L FFT bins, where L is given by: 



L = 




(3.2) 



where: 

f 0 = the nominal transmitter frequency (center frequency); 
B = the total processed bandwidth. 
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2 . 



Derivation of PRF constraints 



To get a moving target in clutter free zone after DFT processing, the motion induced 
doppler shift should be such that shift L exceeds the clutter extent, 

L > NxAf. (3.3) 



From Equations 3. 1, 3.2 and 3.3 



2 NTv 

X 



> NxAf. 



From above v_. in is obtained as 




IT T 2NAR 
X cx 1 

2T2N~AR' 



(3.4) 



(3.5) 



The constraint for maximum velocity can be derived by considering the worst case 
scenario for a target located at the end of the range gate (of width ct/2) at the maximum 
speed. To keep the doppler unambiguous, doppler shift and spread together should not 
exceed the PRF and therefore 

L+P < N-NxAf. (3.6) 

Substituting the value of L and P into this last equation. 
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2NTv vNT ,,,, Ar . 
+ < N( 1-xA/) 

X AR 

2 Tv vT . >r 
+ < 1-tA/ 

X AR 

vT X . X 1 

v + < — -xAf — . 

AR2T 2T 2 T 



(3.7) 



Last term in right hand side of above equation can be recognized as the minimum velocity 
from Equation 3.5 and therefore 

(, X \ X 
( 2AR) 2 T mm 

(3.8) 

X _ X , 

"2 T 2 Jr ' 



1 + - 



2 AT?, 



+ v 



Now, solving for f;. 



2 

> — 


fv I 


f, X \ 


X 




[ 2AR j 



(3.9) 



Equality will give the minimum PRF ((,J. If the factor A/(2AR) «1, then the previous 
equation can be reduced to 




(3.10) 



3. Design Method 

The initial specifications are the nominal carrier frequency (f 0 ), range resolution (A/?), 
maximum radial velocity ( v max ), minimum radial velocity (v mm ), time on target (tot) and due 
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to processing constraints, the maximum number of pulses (N max ). The waveform parameters 
that we want to calculate are the PRP (/ r ), number of pulses (AO, frequency step size (4/0 
and the pulse width (r). 

The proposed method consists of five steps, as follows: 

1. Selection ofPRF from given equation . This equation gives the minimum PRF and 
the exact value can be determined using an iterative procedure later. 

2. The number of pulses is determined from tot and f r : 



N i (tot)f r 





\ 




fr > 



where d B is the antenna beamwidth and d s is the scan rate. 

3. The frequency step size is calculated from 



(3.11) 



A/ = 



c 

2NAR ’ 



(3.12) 



where c is the speed of light (3 x 10 8 m/sec). 

4. The pulse width is chosen from Equation 3.5, which can be rewritten as 



x < .. 

* /, A/ 



(3 13) 



N and r will impact the S/N ratio which determines the radar performance, that is 
probability of detection and probability of false alarm for given radar parameters. Therefore 
it should be ascertained from the following equation that chosen N and r yield adequate S/N 
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ratio to satisfy the radar performance requirements. 



5 P t G 2 X 2 o(xN) 

N (4n) 3 R 4 kT 0 FL' 



where 

P, = the transmitted peak power (watts); 

G = the antenna gain of the transmitter/receiver; 

A = the wavelength (meters); 

< 7 = the radar cross section of the target (square meters); 
r= the pulse width (meters); 

N= the number of pulses coherently integrated within one scan; 

R = the detection range of the target (meters); 
k = Boltzmann’s constant (1.38 x 10' 23 watt-second/°K); 

T 0 = the noise temperature (°K ); 

F = the receiver noise figure; 

L- a loss factor incorporating all system losses. 

By convention, T 0 is taken to be 290°K, which is close to room temperature (300°K) 
and conveniently makes the product kT 0 a round number (4 x 10 ' 21 watt-second/°K). 
5. Flexibility in PRF constraint can be used to satisfy other design objectives such as 
the search range; 
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Search range = number of R u 's . R u 




(3.15) 



k c 



W 2A/' 



where ff is the unambiguous range and B r is the intermediate frequency bandwidth. 
N,f r , rand Af are changed iteratively such that the required search range is obtained. 
Extent of If should not increase clutter as range resolution is already fixed. However 
the number of R u ’s will increase clutter (which should not matter if cancellation is 
adequate). 

The initial specifications and the final calculated design parameters can be summarized 
in the following table: 



Specifications 



Calculated Design Parameters 



Range resolution (AR) 



Minimum PRF (f r ) 



Nominal carrier frequency (f 0 ) 



Number of pulses (N) 



Time-on-target (tot) 



Frequency step size (Af) 



Minimum radial velocity (v mm ) 



Pulse width ( v) 



Maximum radial velocity (v max ) 



Maximum number of pulses (N max ) 



Table 1. Summary of specifications and design parameters. 
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4. Graphical implementation of the Design Method 

Using the equations mentioned in the previous section, a graphical method to design 
the waveform was developed. The first graph is a plot of the range resolution versus minimum 
PRF, for four different nominal carrier frequencies. 




Figure 5. Minimum PRF versus range resolution for four different nominal carrier 
frequencies. 

This graph was created using Equation 3.9 and it can be observed that as the value 
of the desired range resolution decreases (corresponding to an improvement in the range 
resolution), the minimum PRF increases. This graph is used to determine the PRF for a given 
range resolution and carrier frequency. 
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The next graph. Figure 6, is a plot of the PRF versus the number of pulses, and was 



generated using Equation 3.11. 




PRF (KHz) 

Figure 6. Number of pulses versus PRF for three different times-on-target (5, 10, 
and 1 5 ms). 



This plot is used to determine the number of pulses for a given time-on-target (tot) 
and PRF, computed in the previous step from Figure 5. Then, and using Equation 3.12, the 
frequency step size was plotted versus the number of pulses, for three different values of 
range resolution (0.3, 0.5, and 1 meter) as can be observed in Figure 7. This graph is used 
to determine the frequency step for a given range resolution and number of pulses computed 
in the previous step from Figure 6. The frequency step size required for the same number of 
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pulses increases as the desired range resolution decreases. Therefore, the better range 
resolution is desired, the higher the frequency step size must be. 




Figure 7. Frequency step versus number of pulses for three different range resolutions 
(0.3, 0.5, and 1 meter). 



Finally, the last step is to compute the pulse width using Equation 3.13, and all the calculated 
design parameters from Table 1 will be obtained. 

Note that in all the calculations and plots above described v mm was considered to be 
1 50 m/s and v max , 1 000 m/s. The value considered for N max was 2048. 

Consider an example to determine the design parameters for a step frequency radar 
with a nominal carrier frequency of 10 GHz, time-on-target of 15 ms and a desired range 
resolution of 0.3 m. From Figure 5 a minimum PRF of 80 KHz would be obtained. Using this 
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value and Figure 6, N would be approximately equal to 1200. Now with N=1200, dr=0.3m 
and using Figure 7, the frequency step would be approximately 0.42 MHz. The last step is to 
compute the value of the pulse width using Equation 3.14 and the result would came up to 
be t=298 ns. 

Interpolation should be used in this plots whenever the specifications are different 
from the values in the plots. However some times the desired values are beyond the plotted 
values. Additionally, not only we want to compute the desired parameters that satisfy the 
initial specifications, but also we may want to satisfy some parameter constraints. Therefore 
a MATLAB code was written in order to fulfill this requirement. This program can be found 
in Appendix B. 

5. Computer implementation of the Design Method 

The waveform design parameters are constrained to the following values: 

- Pulse width - between 50 nsec and 5000 nsec 

- Frequency step - should be 0.4, 0.6, 0.8 or 1 MHz 

- PRF - should be an integer multiple of 5 KHz 

The program is written in a way so that it will ask for the desired range resolution, nominal 
carrier frequency, and the time-on-target and these data should be entered from the keyboard. 
With this input data the program will compute the minimum values using the equations 
previously presented. These values are iterated until the constraints mentioned above are all 
satisfied and the desired design parameters obtained. The minimum values are shown in 
brackets. The program will also compute the actual range resolution and time-on-target 
corresponding to the rounded values of the PRF, frequency step, pulse width and number of 
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pulses. Figure 8 illustrates how the screen will appear. 



Enter the desired range resolution in meters : 0.3 


Enter the nominal carrier frequency in GHz 


: : 10 


Enter the time-on-target in msec : 15 




The calculated design parameters are : 


(minimum values) 


-PRF = 140.00 KHz 


(80.00 KHz) 


- Frequency step = 0.40 MHz 


(0.42 MHz) 


- Pulse width = 180.00 nsec 


(300.00 nsec) 


- Number of pulses = 2048 


(1200) 


- Actual range resolution = 0.19 m 




- Actual time-on-target = 14.63 msec 





Figure 8. Example of the screen when running program 2. m . 
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IV. STEP FREQUENCY WAVEFORM ANALYSIS USING THE 
AMBIGUITY FUNCTION 

A. INTRODUCTION 

In this chapter, the ambiguity function for a specific step frequency waveform will be 
determined. The waveform with the following parameters is of practical interest: 

t s = 0.1 psec, 

T = 5 psec, 

Af= 1 MHz, 

N = 500, 

where t s is the pulse width, T is the PRI, Af is the frequency step, and N is the number of 
pulses. However, it is difficult to compute and plot this ambiguity function with the desired 
resolution because of the very large amount of computations involved. Not only the amount 
of computations is large but also the plotting is hard to handle because the files can easily 
reach several megabytes. Therefore, waveforms with parameters similar to the waveform of 
interest but easier to compute will be investigated, to bring out the key characteristics of the 
waveform under investigation. First, theoretical dimensions will be verified with a few cases. 
Then, these theoretical dimensions will be used for a particular case of interest which is 
difficult to compute directly. 

Since the ambiguity function of the step frequency waveform contains elements of 
ambiguity functions of linear frequency modulated (LFM) pulses and train of constant 
frequency pulses, the ambiguity functions of these waveforms will be introduced first, along 
with the general introduction to ambiguity function. 
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B. DEFINITION OF AMBIGUITY FUNCTION AND ITS PROPERTIES 
1. Definition of Ambiguity Function 

A radar waveform’s ambiguity function is probably the most complete statement of 
the waveform’s inherent performance. It is a formula that quantitatively describes the 
interference caused by a point target return located at a different range and velocity from a 
reference target of interest. It reveals the range-doppler position of ambiguous responses and 
defines the range and doppler resolution. [Ref. 2, pp. 74] 

This quantitative description of the ability of a waveform to resolve two or more radar 
reflectors at arbitrarily different ranges and velocities, constitutes an important feature for 
quick assessment of the interference level with which a target of interest must compete when 
it is in the vicinity of the other radar reflectors. Although it is seldom used as a basis for 
practical radar system design, it provides an indication of the limitations and utility of 
particular classes of radar waveforms, and gives the radar designer general guidelines for the 
selection of suitable waveforms for various applications. 

The ambiguity function of the waveform s(t) can be defined in terms of the cross- 
correlation of a doppler- shifted version of the waveform, that is s(t) exp (j2^f d t) with the 
unshifted waveform. Using the definition of cross-correlation, it follows that 



where z is the delay time and f d is the doppler frequency shift. Rearranging the terms in the 
integral produces a common form of the ambiguity function as shown in Equation 4.2. 




(4.1) 
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-t-oo 



(4.2) 



= If 5(0 e’ 2 " f ‘'dt\ . 



A normalized expression is obtained by requiring that 



J|s(0|V(=i . 



(4.3) 



With this normalization, the magnitude of the ambiguity function has a value of unity at the 
origin. 

2. Properties of the Ambiguity Function 

The ambiguity function has the following properties: 



Peak value of \yfa 7 f d )\ = |x(0,0 )| = E , (4.4) 

IX(-V/„)I = \X( X ,L )\ , (4.5) 

-t-oo 

lx (T, °)l = 1/ 5 (0 5 *(*-*) dt I > (4.6) 

— oo 
-hoo 

lx(°,/rf)l = I f s 2 (0 e l2%U dt\ , ( 4 . 7 ) 

— OO 

-t-oo -t-oo 

f flx^fjl^df, = E . (4,8) 

— oo — OO 



Equation 4.4, states that the peak value of the ambiguity function occurs at the origin and it 
is equal to E, the energy contained in the echo signal. Equation 4.5 shows the ambiguity 
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function’s symmetry. Equation 4.6 indicates that the ambiguity function along the time delay 
axis is the auto-correlation function of the complex envelope of the transmitted signal. 
Equation 4.7 states that along the frequency shift axis the ambiguity function is proportional 
to the spectrum of s 2 (t). Finally, Equation 4.8 states that the total volume under the ambiguity 
surface is a constant equal to E. 

3. Ambiguity Diagram 

It is common to refer to as the ambiguity surface of the waveform. The 

shape of this ambiguity surface depends entirely on the waveform parameters. The plot of the 
ambiguity surface is called ambiguity diagram. The ideal ambiguity diagram consists of a 
single spike of infinitesimal thickness at the origin and is zero everywhere else. The single 
center spike eliminates any ambiguities, and its infinitesimal thickness at the origin permits the 
frequency and the echo delay time to be determined simultaneously to an high degree of 
accuracy. It also permits the resolution of two targets no matter how close together they are 
on the ambiguity diagram. Naturally, this ideal diagram does not exist. The two reasons for 
this can be found in the properties: first, and accordingly with Equation 4.4, the maximum 
height of the ambiguity function is E and secondly the volume under the surface must be finite 
and equal to E, as stated by Equation 4.8. However, a reasonable approximation is given in 
Figure 9. This ambiguity function only has one peak and therefore does not cause any 
ambiguity. However the single peak might not be narrow enough to satisfy the requirements 
of accuracy and resolution. If the single central peak is made too narrow, it may cause other 
smaller peaks to occur in regions other than the origin, and therefore cause ambiguities. The 
requirements for accuracy and unambiguity are not always possible to satisfy simultaneously. 
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Figure 9. An approximation to the ideal ambiguity diagram “After Ref. [1].” 

The particular waveform transmitted by a radar is chosen to satisfy the requirements 
for (1) detection, (2) measurement accuracy, (3) resolution, (4) ambiguity, and (5) clutter 
rejection. The ambiguity diagram may be used to assess how well a waveform can achieve 
these requirements. Each of these will be discussed briefly. 

The requirements for detection do not place any demands on the shape of the 
transmitted waveform except that it be possible to achieve with practical radar transmitters, 
and the maximum value of the ambiguity function is an indication of the detection capabilities 
of the radar. The accuracy with which the range and the velocity can be measured by a 
particular waveform depends on the width of the central spike along the time and frequency 
axis. The resolution is also related to the width of the central spike, but in order to resolve 
two closely spaced targets the central spike must be isolated. It cannot have any high peaks 
nearby that can mask another target close to the desired target. A waveform that yields good 
resolution will also yield good accuracy, but the reverse is not always so. 
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A continuous waveform (a single pulse) produces an ambiguity diagram with a single 
peak. A discontinuous waveform can result in peaks in the ambiguity diagram at other values 
of x,f d . The pulse train is a common example. The presence of additional spikes can lead to 
ambiguity in the measurement of target parameters. An ambiguous measurement is one in 
which there are several choices available for the correct value of a parameter, but only one 
choice is appropriate. Thus the correct value is uncertain. The ambiguity diagram permits a 
visual indication of the ambiguities possible with a particular waveform. The ambiguity 
problem, detection and accuracy are related to a single target, whereas resolution applies to 
multiple targets. 

The ambiguity diagram may be used to determine the ability of a waveform to reject 
clutter by superimposing on the t, f d plane the regions where clutter is found. If the 
transmitted waveform is to have good clutter-rejection properties the ambiguity function 
should have little or no response in the regions of clutter. 

The problem of synthesizing optimum waveforms based on a desired ambiguity 
diagram specified by operational requirements is not normally feasible. The approach to 
selecting a waveform with a suitable ambiguity diagram is not systematic but rather by trial 
and error. 

The name ambiguity function is somewhat misleading since this function describes 
more about the waveform than just its ambiguity properties. This name was given to this 
function in order to demonstrate that the total volume under it, is a constant equal to E, 
independent of the shape of the transmitted waveform. Thus the total area of ambiguity, or 
uncertainty, is the same no matter how the ambiguity surface is distributed over the x ,f d 
plane. [Ref. 1, pp. 418-420] 
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C. AMBIGUITY FUNCTION OF A SINGLE PULSE 



The unmodulated single pulse is widely used in old generation radars for search and 
track functions (as magnetrons are not easily modulated) and where range accuracy and 
resolution requirements can be met with a pulse wide enough to provide sufficient energy for 
detection. It has the minimum ratio of time sidelobe extent to compressed pulse width and is 
used in inexpensive radars where signal generation and processing costs must be minimized. 

The single pulse of a sine wave can be defined as 

S(t)=s(t)e J2ltfo ‘ , (4.9) 



where s(t) is the complex envelope of the signal, defined as follows: 

s(t) = 1, if 0< t< t s , 

= 0, elsewhere, (4.10) 

with t s equal to the pulse width. 

Using Equation 4.2, the ambiguity function of the single pulse can be written as 
follows: 






(i-Ml 


sin[7t^(l-|T|/g/J 


( 'J 


-W 5 )f d 



zero elsewhere. (4. 1 1) 



The figures for this section are in Appendix A. The ambiguity diagram of the single 
pulse is shown in Figure A. 1 and the contour plot can be observed in Figure A.2. Profiles of 
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the ambiguity function, taken through its peak, are shown in Figures A.3 and A.4. The time 
profile is a triangle (the autocorrelation of a rectangle is a triangle) with a half-voltage width 
equal to the pulse width. The frequency profile is a sine function whose main lobe width is 
twice the reciprocal of the pulse width. 

D. AMBIGUITY FUNCTION OF THE LINEAR FREQUENCY MODULATED 

PULSE 

The linear frequency modulated (LFM) pulse is commonly used to increase range 
accuracy and resolution when long pulses are required to get reasonable signal-to-noise ratios 
(10 to 20 dB). This waveform can be used for detection of targets with unknown velocity 
since the doppler sensitivity is low. 

The LFM pulse can be represented mathematically as 



where s(t) is the same as defined in Equation 4 . 10 , and k is the rate of frequency change in 
Hz/sec. The ambiguity function of the LFM pulse can be written as in Equation 4.13. 




(4.12) 




- \fso) sxt-^^Wd^ 



sin [71^(1 -|'c]/( s )(/ t/ + £'c)] 




zero elsewhere. (4.13) 
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Comparing this equation with Equation 4.2, we can see that they are identical, except that f d 
is replaced by f d - kt. Therefore , we can conclude that the ambiguity function of the LFM 
pulse is a shifted version of the ambiguity function of a single pulse. This relation is 
represented graphically in Figure 10 and it can be seen that the ambiguity function of the LFM 
pulse is just a rotated version of the ambiguity function of the single pulse. 



Figure 10. Contour comparison for the single pulse and the LFM pulse. 

The figures for this section are in Appendix A. Figures A. 5, A. 6, A. 7, and A. 8, 
represent the ambiguity diagram, contour plot, and time and frequency profiles (taken through 
the peak), respectively. 
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E. AMBIGUITY FUNCTION OF THE CONSTANT FREQUENCY PULSE TRAIN 



Practical radars employ waveforms consisting of constant frequency pulse trains and 
therefore it is important to study their ambiguity functions. The ambiguity functions for these 
pulse trains reveal ambiguous responses in range and doppler. Performance for specific 
surveillance applications can be understood in terms of unambiguous range-doppler regions 
of operation determined by radar pulse repetition frequency, pulse duration, and pulse 
bandwidth. For high resolution applications, the ambiguity surface of individual pulses of the 
pulse train is also of interest. [Ref. 2, pp. 74-75] 

For single pulse delay and frequency measurement accuracies depend on the single 
parameter of pulse width. With the pulse train this situation can be avoided. The delay 
accuracy depends on the pulse width as before, but the frequency accuracy is now determined 
by the total duration of the pulse train. Therefore, both accuracies are independent of one 
another. The price that has to be paid for this capacity of independently controlling delay and 
frequency accuracies is that additional peaks occur in the diagram, which in turn cause range 
and doppler ambiguities. In practice, the radar designer tries to select the PRI in order to 
make all targets of interest appear in the vicinity of the central peak, and all the other peaks 
occur as far from this region as possible. 

Most radars use this type of waveform, which can be mathematically represented as 
follows: 



AM 



S(i) = 



J2*f 0 t 



n=0 



(4.14) 



where s(t) is the complex envelope of the single pulse of the transmitted signal, defined as 
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s(t-nT) = 1, 
= 0 , 



if nT< t< nT+t s , 



elsewhere. 



(4.15) 



In this equation, tV is the number of pulses, T is the PRI and t s is the pulse width. The 
ambiguity function of the constant frequency pulse train can be written as 



Figure 1 1 is a level contour of the ambiguity surface of a constant frequency pulse 
train. This figure shows that the width of the central peak along the delay axis is 
approximately twice the pulse width, and that along the frequency axis is 2/NT. The 
interpeak distance along the delay axis is equal to the PRI, and along the frequency axis 
1/PRI, which is the PRF of the waveform. The total non-zero extent along the delay axis is 
2NT and along the frequency axis is 2/t s . Figure 12 represents an example of an ambiguity 
diagram for the constant frequency pulse train. The contour plots in Figures 13 and 14 show 
that the results obtained agree with the theoretical ones previously defined. This can also be 
observed in Figures 15 and 16, which represent the time and frequency profiles, respectively. 
The time profile of the central pulse is a triangle and the frequency profile is a sine function, 
as expected. 




(4.16) 



Changing variables ( t-nT=t ' ), Equation 4. 16 becomes 



TV- 1 TV- 1 

\xt,rn - lE E ‘ mr ‘" T 




39 



f d 





Figure 1 1 . Level contour of the ambiguity surface of a constant frequency pulse train. 
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Figure 12 Ambiguity diagram of a constant frequency pulse train 
(N~5, pulse width=l microsec, PRl=5 microsec, Af=0) 
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Figure 13. Contour plots of the ambiguity diagram of a constant frequency pulse train 
(N=5, pulse width=l microsec, PRI=5 microsec, Af=0) 

(a) Global view. 

(b) Magnified view. 
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frequency [MHz] frequency [MHz] 





Figure 14. Contour plots of the ambiguity diagram of a constant frequency pulse train. 
(N=5, pulse width =1 microsec, PRI=5 microsec, Af=0) 

(a) Magnified view. 

(b) Magnified view of the central peak. 
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Figure 15. Time profiles of the ambiguity diagram of a constant frequency pulse train. 
(N=5, pulse width =1 microsec, PRI=5 microsec, Af=0) 

(a) Global view. 

(b) Magnified view. 
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(b) 



Figure 16. Frequency profile of the ambiguity diagram of a constant frequency pulse 
train . (N=5, pulse width =1 microsec, PR1=5 microsec, Af=0) 

(a) Global view. 

(b) Magnified view. 
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F. AMBIGUITY FUNCTION OF THE STEP FREQUENCY WAVEFORM 

The transmitted signal of the step frequency radar (as shown in Figure 1) can be 
represented mathematically as follows: 

S t (t ) = A t £ s(t-nT) e J2n ^ nAf)t , (4. 1 8) 

n-0 



with s(t) defined in Equation 4.15, and n = 0 , ... , N-l. A, is the amplitude of the transmitted 
signal. Equation 4. 18 can be rewritten as follows: 



S t (t) = A, 



N - 1 

^2 s(t-nT) e J2nn ^ 1 

n=0 



J 2 *fo‘ 



= A t 5 (0 e 



j 2it/o« 



(4.19) 



The expression in brackets is the complex envelope of the pulse sequence represented by S(t). 
The received signal can be written as 

S r (t ) = A r S(t-x) e j2n(foVdK, ~ x} . (4.20) 



Substituting the complex envelope S(t) into Equation 4.2, will give an expression for 
the ambiguity function 



Xi'tJd) I = I f 



N-\ 



AM 



J2s(t-mT) Y2 s*(t-nT-x) 

m~ 0 n =0 



e -j2nnLf(t-x) 



e 2nfd ’dt | . (4.21) 
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Changing variables ( t-mT = t ’ ), the final expression for the ambiguity function of the step 
frequency waveform is obtained: 

N - 1 N - 1 

g)2v.m 2 A/T gjlKmfjT ~j2nmnA/T ^j2nnAfz 

m-0 «=0 

(4.22) 

0 s\t'-{n-m)T-z) e^m-'OA/r' e iW dt' I . 

The theoretical dimensions of the contours of the ambiguity surface for the step 
frequency waveform are shown in Figure 17. The overall dimensions are of length 2NT along 
the delay axis and 2NA/ along the frequency axis. It can be observed that the distance 
between component contours along the delay axis is equal to the PRI of the waveform and 
along the frequency axis is equal to the inverse of the PRI. The ambiguity diagram of the step 
frequency waveform can be obtained from the ambiguity diagram of the constant frequency 
pulse train (same parameters), by rotating the delay axis by A/7T. This can be stated more 
accurately by saying that the ambiguity function of the step frequency radar is x( T , (4/7T) f d ) 
where x(x, f d ) is the ambiguity function of the constant frequency pulse train. Apart from 
rotation of horizontal axis (and all cuts parallel to it) the individual spikes also rotate. This is 
similar to linear frequency modulation pulses. The extent of the central column of spikes along 
the frequency axis for zero delay is 2/t s , value which decreases when we move away from the 
center. The central spike as shown in Figure 17 is inclined A/7T with respect to the horizontal 
axis, indicating the range-doppler coupling of the step frequency waveform. The projection 
of the central peak on the delay axis is 2t s , which is the overall delay uncertainty. However 
if the doppler is known, the delay uncertainty is given by 2/(NA f) as indicated in the figure. 
The projection of the central peak on the frequency axis is also indicated in the figure. If the 




WJ d )\ = iE E 
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Figure 17. Level contour of the ambiguity surface of a step frequency waveform. 
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delay is known, the frequency uncertainty is 2/NT which is the same as for the constant 
frequency pulse train. However if target delay is not known, frequency uncertainty is much 
larger due to range-doppler coupling. 

G. AMBIGUITY FUNCTION OF STEP FREQUENCY WAVEFORMS 

Four waveforms will be studied in this section. These waveforms are defined by the 
following parameters: 

1. N=5, t s =1 ps, PRI =5 ps, Af =1 MHz. 

2. N=10, t s =1 ps, PRI =5 ps, Af =1 MHz. 

3. N=10, t s =0.1 ps, PRI =5 ps, Af =1 MHz. 

4. N=500, t s =0.1 ps, PRI =5 ps, Af =1 MHz. 

1. First waveform 

The results for this waveform are shown in Figures 18 through 22. The ambiguity 
diagram is represented in Figure 18. It can be observed that it is spiky like the one for the 
constant frequency pulse train, but now it is tilted at an angle. The global view of the contour 
plot of the ambiguity function for this waveform is shown in Figure 1 9(a). The rotation of the 
plot with respect to the axes is obvious in the figure and the rotation angle is given by A/7T. 
The delay and frequency axis dimensions match the theoretical values as from -NT to NT and 
from -NA/ to NA f . Figure 19(b) is a close up view of the contour plot. Figure 20(a) is a 
further magnification of the contour plot and Figure 20(b) shows the details of the central 
peak. From the contour plot it is clear that the central peak is a inclined ridge with the highest 
magnitude in the center of the contour. Figure 21(a) is a cut of the ambiguity surface along 
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Figure 1 8 Ambiguity diagram of a step frequency waveform 

(N=5, pulse width=l microsec, PR1=5 microsec, Af=lMHz) 
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Figure 19. Contour plots of the ambiguity diagram of a step frequency waveform. 
(N=5, pulse width=l microsec, PRI=5 microsec, Af=l MHz) 

(a) Global view. 

(b) Magnified view. 
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Figure 20. Contour plots of the ambiguity diagram of a step frequency waveform. 
(N=5, pulse width =1 microsec, PRI=5microsec, Af=l MHz) 

(a) Magnified view. 

(b) Magnified view of the central peak. 
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Figure 21. Time profiles of the ambiguity diagram of a step frequency waveform. 
(N=5, pulse width =1 microsec, PRI=5 microsec, Af=lMHz) 

(a) Global view. 

(b) Magnified view. 
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Figure 22. Frequency profile of the ambiguity diagram of a step frequency waveform. 
(N=5, pulse width =1 microsec, PRi=5 microsec, Af=lMHz) 

(a) Global view. 

(b) Magnified view. 
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the delay domain at f d - 0, which portrays only one peak. This suggests that there is no range 
ambiguity for the waveform if the target doppler is known. Figure 21(b) is a magnification of 
Figure 21(a) which gives null-to-null width of 0.4 which is the same as computed from the 
formula 2/(NA/). Note that this is also the output of the matched filter in the time domain. 
The range resolution of this waveform is 0.2 as given by 1/(NA/). Figure 22(a) is a cut of the 
ambiguity surface along the frequency domain at t= 0. Theory suggests that there should be 
(2T/t s )-l peaks in the frequency domain, which is matched by the actual value of 9. Figure 
22(b) is a detail of the central peak and should have a null-to-null width of 2/(NT) which 
matches with the actual value of 0.08. 

2. Second waveform 

For the second case, the waveform has the same parameters as in the previous one 
except that the number of pulses is increased to 10. This would affect the range resolution 
(1/NAf) and the doppler resolution (1/NT) which are now 0.1 and 0.02 respectively. 
However, some other quantities of interest such as the component contours in the frequency 
axis stay the same. The results can be observed in Figures 23 through 26. Figure 23(a) gives 
the global view of the contour plot of the ambiguity diagram for this waveform. Figures 
23(b), 24(a) and 24(b) are successively increasing magnifications of the contour plot. Figure 
24(b) is a detailed plot of the central peak giving a range resolution and frequency resolution 
of 0. 1 and 0.02, respectively. This can be confirmed in Figures 25(b) and 26(b) which are 
magnified views of the time and frequency profiles, respectively. 
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Figure 23. Contour plots of the ambiguity diagram of a step frequency waveform. 
(N=10, pulse width=l microsec, PRI=5 microsec, Af=l MHz) 

(a) Global view. 

(b) Magnified view. 
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Figure 24. Contour plots of the ambiguity diagram of a step frequency waveform. 
(N=10, pulse width=l microsec, PRI-5 microsec, Af=l MHz) 

(a) Magnified view. 

(b) Magnified view of the central peak. 
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Figure 25. Time profiles (cut at zero frequency) of the ambiguity diagram of a step 
frequency waveform. (N=10, pulse width=l microsec, PRI=5 microsec, 
Af=lMHz) 

(a) Global view. 

(b) Magnified view. 
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Figure 26. Frequency profiles of the ambiguity diagram of a step frequency waveform. 
(N=10, pulse width =1 microsec, PRI=5 microsec, Af=lMHz) 

(a) Global view. 

(b) Magnified view. 
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3. 



Third waveform 



For this waveform the pulse width is reduced to 0. 1 ps and the rest of the parameters 
are the same as in the previous case. The results can be observed in Figures 27 through 32. 
Figure 29 represents the central peak of the ambiguity diagram of this step frequency 
waveform which is similar to the one in Figure A. 5. 

Now that all the theoretical dimensions previously defined were verified, they can be 
used to calculate the dimensions of any particular case of interest. This will be done for the 
fourth and last case. 

4. Fourth waveform 

Using the results from Section F, which have been confirmed for the waveforms 
discussed so far, ambiguity diagram figures were sketched for a waveform with 500 pulses. 
The amount of computations required to generate the ambiguity surface by computer is 
enormous. The 3D ambiguity diagram of this waveform is several magnitudes more complex 
than 3D pictures of midtown Manhattan. One can realize the problem of resolution of this 
figure and its printing on a normal size paper or even in a much larger size of paper. 

A contour plot of the ambiguity diagram is sketched in Figure 33. Each short slant line 
is a spike. In doppler dimension the spikes are PRF (200 KHz) apart and spread over 2/t s (20 
MHz). In delay dimension (at f d =0) there are about ten significant spikes (in reality there are 
more but others will be small in magnitude). This means that the IF filter will pass ten 
frequency lines, which implies that the return from ten different ranges due to separate pulses 
may arrive at the same time. Thus, one may say that the target will be ambiguous in range and 
that it can lie in any of the ten range zones. Apart from the fact that the target range ambiguity 
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Figure 27. Contour plot of the ambiguity diagram of a step frequency waveform. 
(N=10, pulse width=0.1 microsec, PRI=5 microsec, Af=l MHz) 




delay [microsec] 



Figure 28. Time profile (cut at zero frequency) of the ambiguity diagram of a step 
frequency waveform. (N=10, pulse width=0.1 microsec, PRI=5 microsec 
Af=l MHz) 
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Figure 29. 3D plot of the central peak of the ambiguity diagram of a step frequency 
waveform. (N=10, pulse width=0.1 microsec, PRI=5 microsec, Af=l MHz) 




delay [microsec] 



Figure 30 Contour plot of the central peak of the ambiguity diagram of a step 
frequency waveform. (N=10, pulse width=0.1 microsec, PRI=5 microsec, 
Af=l MHz) 
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Figure 3 1 . Time profile (cut at zero delay) of the central peak of the step frequency 
waveform. (N=10, pulse width = 0.1 microsec, PR1 = 5 microsec, 
Af=lMHz) 




frequency [MHz] 



Figure 32. Frequency profile (cut at zero delay) of the central peak of a step 
frequency waveform. (N=10, pulse width=0.1 microsec, PRI=5 microsec, 
Af=l MHz) 
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Figure 33. Contour plot of the ambiguity diagram of a step frequency waveform. 
(N=500, pulse width=0. 1 microsec, PRI=5 microsec, Af=l MHz) 
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needs to be resolved it will also compete against clutter from ten different range zones. 

If the waveform parameters are picked up judiciously, the doppler from a target 
moving at the highest expected velocity may be kept unambiguous. In such cases, spikes 
along the doppler dimension may be ignored, otherwise they have to be considered for clutter 
calculations. Figures 34 and 35 represent cuts along the delay and frequency dimensions 
indicating the potential range and doppler ambiguities, respectively. Figure 36 gives the detail 
of the central peak (note that other peaks or spikes are of similar shape but may be of 
different magnitude and located at different places in range-doppler map). At f d =0 the peak 
is 4 ns wide (null-to-null) in delay dimension. Similarly, the width of the peak in frequency 
dimension (at t= 0) is given by 2/NT which comes out to be 800 Hz. 
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Figure 34. Time profile (cut at zero frequency) of the ambiguity diagram of a step frequency 
waveform. (N=500, pulse width=0.1 microsec, PRI=5 microsec, Af=l MHz) 
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Figure 35. Frequency profile (cut at zero delay) of the ambiguity diagram of a step frequency 
waveform. (N=500, pulse width=0. 1 microsec, PR1=5 microsec, Af=l MHz) 
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Figure 36. Dimensions of the central peak for the case N=500, pulse width = 0.1 (is, 
PRI = 5 us, Af = 1 MHz. (Values indicated are null-to-null; generally 3dB 
values are given which are a little less than half of the indicated.) 
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V. CONCLUSIONS 



This thesis investigates the use of the step frequency waveform, its design and 
analysis using the ambiguity function. The step frequency waveform achieves a high range 
resolution by coherently processing the returns from N transmitted pulses, each with a 
different carrier frequency, monotonically increased from pulse to pulse. This waveform has 
two main advantages over the conventional high resolution waveforms. First, the narrow 
instantaneous bandwidth does not require a high A/D sampling rate which can be a major 
limitation in the system design and secondly, the step frequency waveform can be 
implemented on existing radar equipment using a step frequency synthesizer. 

A design procedure for detection of small targets with a surface (land or sea) based 
step frequency radar which employs a high PRF waveform was developed. The proposed 
method determines the waveform parameters for given radar specifications. Two different 
implementations of this method were presented: the graphical and the computer 
implementation. The graphical implementation is a simple and quick way of calculating the 
minimum design parameters. The computer implementation considers the parameter 
constraints that are present when designing a waveform and associated equipment. 

The ambiguity function, an important tool for the analysis of any waveform, was 

defined along with its properties. A practical step frequency waveform may include hundreds 

of pulses. It is difficult to compute the ambiguity function of a waveform with a large number 

of pulses. Waveforms with parameters close to the waveform of practical interest but easier 

to compute were investigated, in order to bring out the key characteristics of the waveform 

of interest. The theoretical dimensions of the step frequency waveform were defined and 
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verified with a few different cases. Once this was done, any particular case of interest could 
be investigated. The ambiguity diagram of one step frequency waveform of interest was 
sketched. The ambiguity function of the step frequency waveform contains elements of the 
ambiguity function of a single LFM pulse and train of constant frequency pulses. The 3D plot 
of the ambiguity function of a step frequency waveform is spiky like the one of a constant 
frequency pulse train and is tilted at an angle like LFM. The width of the central spike along 
the delay axis at zero frequency is 2/(NA$ as compared to 2t s for the constant frequency 
pulse train, therefore allowing to decrease the width of the spike by increasing the product 
NA/. This permits increasing the range resolution without increasing the instantaneous 
bandwidth. For the constant frequency pulse train, high range resolution can only be achieved 
if the pulse width is decreased. This will cause the increase of the instantaneous bandwidth 
and therefore a higher A/D sampling rate. On the frequency axis, the spike width at zero delay 
is equal to 2/(NT) as it also for the constant frequency pulse train case. Therefore, the 
frequency resolution is not improved with the step frequency waveform. 

Further work on the comparative analysis of step frequency waveform should be 
performed to compare it with conventional high PRF and medium PRF waveforms for 
detection of small targets. Its use in inverse synthetic aperture radar (ISAR) for target 
identification should also be investigated. 
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APPENDIX A 



This appendix contains the figures for Chapter IV, Section C - Ambiguity function of 
a single pulse, which are Figures A. 1 through A.4. It also contains the figures for Section D - 
Ambiguity function of the linear frequency modulated pulse, which are Figures A. 5 through 
A.8. 
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Figure A 1 Ambiguity diagram of a single pulse 
(pulse width= I microsec) 
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Figure A. 2. Contour plot of the ambiguity diagram of a single pulse, 
(pulse width=l microsec) 
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Figure A.3. Time profile of the ambiguity diagram of a single pulse, 
(pulse width =1 microsec) 




Figure A.4. Frequency profile of the ambiguity diagram of a single pulse. 

(pulse width =1 microsec) 
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Figure A 5 Ambiguity diagram of a L.FM pulse 
( pulse width- 1 microsec) 
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Figure A. 6. Contour plot of the ambiguity diagram of a LFM pulse, 
(pulse width=l microsec) 
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Figure A.7. Time profile of the ambiguity diagram of a LFM pulse, 
(pulse width =1 microsec) 




Figure A. 8. Frequency profile of the ambiguity diagram of a LFM pulse. 

(pulse width =1 microsec) 
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APPENDIX B. MATLAB SOURCE CODES 



% File name: programl.m 
% Title: Plot generation 
% Date of last revision : 4 Mar 1 996 
% Written by: Paulo A. Soares 

% Comments: This program generates 3 plots for the SFWF design. 

Oy ^**** ********************************************************* ****** ********* 

clear 

% specifications 

vmin=l 50;vmax=1000; 

fO 1 =3e9;f02=5e9;iD3=7e9;j04= 1 0e9; 

tot 1 =5e-3 ;tot2= 1 Oe-3 ;tot3= 1 5e-3 ;tot4=20e-3 ; 

% computations & plots 

c=3e8; 

lambda 1 =c/fO 1 ;lambda2=c/102 ; 
lambda3=c/fl)3 ;lambda4=c/f04; 

-% plot $1 

dr=0.0 1:0.0 1:1; 

frl=(dr.\vmax)+((2/lambdal )*(vmax+vmin)); 
fr2=(drAvmax)H<(2/lambda2)*(vmax+vmin)); 
fr3=(dr.\vmax)+((2/lambda3)*(vmax+vmin)); 
fr4=(dr.\vmax)+{(2/lambda4)*(vmax+vmin)); 
figure(l) 

plot(dr,ff l/le3,'w’),hold on 
plot(dr,ff2/le3,V),hold on 
plot(dr,fr3/le3,'w'),hold on 
plot(dr,fr4/le3,W) 
axis([0.01, 1,10,1 10]) 
xlabel('Range Resolution (meters)') 
ylabel('Minimum PRF (Khz)') 
text(0.75,27,’fo=3GHz’),text(0.75,42,'fo=5GHz') 
text(0.75,57,’fo=7GHz'),text(0.75,81 ,'fo=l 0GHz') 
text(0.2,100,’fo = nominal carrier frequency') 

% plot #2 

lambda=0.03 :0.00 1 :0. 1 29; 
k=(dr.\vmax); 

fr 1 =(lambda.\2)*(vmax+vmin)+k; 

N1 l=totl*frl; 

N21=tot2*frl; 

N31=tot3*frl; 

figure(2) 

plot(frl/le3,Nl l,'w’)Jbold on 
plot(frl /I e3,N2 1 , V) Jbold on 
plot(fr 1 /le3,N3 1 ,'w') Jbold on 
axis([20, 180,0,3000]) 
xlabel(’PRF (KHz)') 
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ylabel('N') 

text( 1 40,900,'tot=5ms'),text( 1 40, 1 7 00,’ tot= 1 0ms') 
text(l 40,25 50,’tot=l 5ms') 
text(50, 2500, ’tot = time-on-target') 

% plot #3 

deltarl=0.3; 

deltar2=0.5; 

deltar3~l; 

N=l: 1:2048; 

deltafl=(2*N*deltarl * le6).\c; 
delta£2=(2*N*deltar2* le6).\c; 
deltaD=(2*N*deltar3* le6).\c; 
figure(3) 

plot(N,deltafl ,'w'),hold on 
plot(N,deltaf2,'w'),hold on 
plot(N,deltaf3,'w') 
axis([0,2 100,0,1]) 
xlabel('N') 

ylabel(’Frequency Step (MHz)’) 

text(l 700, 0.12, 'dr=lm'),te:rt(l 700,0.21, 'di^.5m’),text(l 700, 0.33, ’dr=0.3m') 
text(l 200,0.8,’dr = range resolution') 
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% Filename: program2.m 
% Title: Main program for SFWF design 
% Date of last revision: 28 Feb 1996 



% Written by: Paulo A. Soares 

% Comments: This program computes the exact design parameters frmin, df, N and tau, and then 
% rounds it in accordance with previously defined criteria. 



clear 

c=3e8;Nmax=2048; 
vmin= 1 50;vmax= 1 000; 

dr=input(' Enter the desired range resolution in meters : '); 
while dr <= 0 

dr=input(' You must enter a positive value. Please try again : ’) 
end 

fl>=input(' Enter the nominal carrier frequency in GHz : f ); 
while fD <= 0 

f0=input(' You must enter a positive value. Please try again : ’) 
end 

lambda=c/(fO* 1 e9); 

tot=input(’ Enter the time-on-target in msec : ’); 
while tot <= 0 

tot=input(’ You must enter a positive value. Please try again : ’) 
end 

frmin=(((2/lambda)*(vmax+vmin))+(vmax/dr))/l e3; 
fr=ceil(frmin); 
while rem(fr,5) ~ 0 
fr=fH-l; 
end 

Nmin=ceil(tot*fr); 

N=Nmin; 
if N<=2 
N=2; 
end 

elseif N=3 | N=4 
N=4; 
end 

elseif N>4 & N<=9 
N=8; 
end 

elseif N>9 & N<=17 
N=16; 
end 

elseif N> 17 & N<=34 
N=32; 
end 

elseif N>34 & N<=68 
N=64; 
end 

elseif N>68 & N<=136 
N=128; 
end 

elseif N> 136 & N<=272 
N=256; 
end 
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elseif N>27 2 & N<=544 
N=512; 
end 

elseif N>544 & N<=1088 
N=1024; 
end 

elseif N>1088 
N=Nmax; 
end 

Nround=N; 
frround=N round/tot; 
finroiind=ceil(frroimd) ; 
while rem(ffround.5) — =0 
fm>und==ffround+ 1 ; 
end 

totround=N round/frround ; 

dfmin^/(2*Nmin*dr* 1 e6) ; 

df=c/(2*Nround*dr* le6); 

stepsizedf=0.2; 

dfround=ceil(df7stepsizedf); 

dffound=dfround*stepsizedf; 

taumin— (2 * viiiiii)/( lambda* dfmin* frmin) ; 

tau=(2*vmin)/(lambda*dffound*fn‘ound); 

stepsizetau=10; 

tauround=ceil(tau/stepsizetau); 

tauround=tauround* stepsizetau; 

acdr=(c/(2*Nround*dfround))* 1 e-6; 

stepsizeacdr=0.0 1 ; 

acdrround=eeil(acdr/stepsizeacdr); 

acdiTound=acdrround* stepsizeacdr; 

fprintf(’\n The calculated design parameters are An’) 

fprintfC (minimum values in brackets)\n\n’) 

fprintfC -PRF = %3.2f KHz (%3 .2f KHz^’.frrounifrmin) ; 

fprintfC -Frequency step = %3.2f MHz (%3.2f MHz)\n^dfround,dfmin); 
fprintfC -Pulse width = %3 .2f nsec (%3 .2f nsec)\n',tauround,taumin); 

fprintfC -Number of pulses = %4.0f (%4.0f)W,Nround,Nmin); 

fprintfC -Actual range resolution = %3 .2f m\n\acdrround); 
fprintfC -Actual time-on-target = %3.2f msec \totround); 
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% Filename: program3.ni 

% Title: Ambiguity function of the simple pulse by correlation method 
% Date of last revision: 6 Apr 1996 
% Written by: Paulo A. Soares 

% Comments: This program draws four plots of the ambiguity function for a single pulse. 
%********************************************************************************* 



clear; 

tau=l; 

nx=50; 

ny=201; 

t=[0:nx-l]*tau/nx; 

dt=t(2)-t(l); 

z=[]; 

fd-linspace(- 1 0/tau, 1 0/tau,ny); 

ul=ones(l,nx); 

for m=l:ny 

u2=u 1 . *exp(j * 2 *pi * fd(m) * t) ; 
c=xcorr(u2,u 1 ). *dt; 
z=[z;abs(c)]; 
end 

t= [flip lr(-t) ,t( 2 :nx)] ; 
figure(l) 
mesh(tfclz),grid 
xlabel('delay [microsec]’) 
ylabel(’frequency [MHz]’) 
zlabel(’magnitude’) 
figure(2) 

contour(tjcLzJ),grid 
xlabel( ’delay [microsec]’) 
ylabel('firequency [MHz]’) 
figure(3) 

plot(U((ny+l)/2,:),’w’) 
xlabel(’delay [microsec]’) 
ylabel(’magnitude’) 
figure(4) 

plottfd^O^nxVw’) 
xlabel(’firequency [MHz]') 
ylabel(’magnitude') 
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% Filename: program4.m 

% Title: Ambiguity function of the LFM pulse by correlation method 
% Date of last revision: 6 Apr 1996 
% Written by: Paulo A. Soares 

% Comments: This program draws four plots of the ambiguity function for a LFM pulse. 

%* *******************************************************************^^*********** 

clear; 

tau=l; 

mu=4; 

nx=50; 

ny=201; 

t=[0:nx-l]*tau/nx; 
dt=t(2)-t(l); 
z =□; 

fd=linspace(- 1 0/tau, 1 0/taiuiy); 
for m=l:ny 

ul=exp(j*pi*mu.*t.*t); 
u2=conj(u 1 ).*exp(j *2*pi*fd(m)*t); 
c=conv(u2iliplr(u 1 )) . *dt; 
z=[z;abs(c)]; 
end 

t=[fliplr(-t),t(2:nx)]; 
figure(l) 
mesh(t/d ? z),grid 
xlabel(’delay [microsec]') 
ylabel(’frequency [MHz]’) 
zlabel(’magnitude') 
figure(2) 

contour( t,fd,z, 5 ) , grid 
xlabel(’delay [microsec]’) 
ylabel(’frequency [MHz]’) 
figure(3) 

plot(fz((ny+l )/2,:)/w’) 
xlabel(’delay [microsec]') 
ylabel('magnitude’) 
figure(4) 

plot(fd,z(:,nx),W) 
xlabel(’frequency [MHz]’) 
ylabel(’magnitude’) 
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% Filename: programS.m 

% Title: Ambiguity function of the step frequency radar by the correlation method 
% Date of last revision: 26 Apr 1996 
% Written by: Paulo A. Soares 

% Comments: This program computes and plots several figures of the ambiguity function of the 
% step frequency radar, using the correlation method. 

%********** ************************ *********************************************** 



%- 

% 

%- 

clear; 

tau=l; 

T=5; 

N=10; 

df=l; 

sampfreq=40; 

NT=N*T; 

nsamptau=sampfreq*tau; 
nsamptO=sampff eq*(T -tau); 
nsamptau=floor(nsamptau); 
nsamptO=floor(nsamptO); 
nsampT=nsamptau+nsamptO; 



Parameters initialization 



% clears all variables 
% tau=pulse width 

% T=period or pulse repetition interval (PRI) 

% N=number of pulses 
% df=frequency step 
% sampfreq=sampling frequency 
% NT=total period of the waveform 
% nsamptau=number of samples in one pulse width 
% nsamptO=number of samples in the interpulse period 



% nsampT=number of samples in one PRI 



%- 

% Definition of axis and resolution 

%- 

P=[ones( 1 > nsamptau) > zeros( 1 ,nsamptO)] ; 

t=linspace(0,NT,N*nsampT); % time axis definition 

ny=30 1 ; % number of points in the freq. axis 

fd=linspace(-5/T,5/T ? ny); % freq. axis definition 

dt=t(2)-t(l); % dt=time interval 

%- 

% Computations 

%- 

s=[|; 

m=l:l:N; 
for n=l :N 

s=[s,exp(j*2*pi*(m(n)-l)*df* : t(nsampT*(n-l)+l:nsampT*n)) *P]; 
end 

z=0; 

for L=1 :ny 
f==fd(L); 

s 1 =conj (s).*exp(j *2*pi*f* c t); 
s2=s; 

c=conv(sl ,fliplr(s2)).*dt; 

%z=[z;(abs(c)). A 2]; 

z=[z;abs(c)]; 

end 

[U]=size(z); 

t=[fliplr(-t) ? t(2 :N*nsampT); 
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Plotting of the results and axis labeling 



%r 

% 

%- 

figure(l) 
mesh(tf(lz),grid 
xlabel('delay [microsec]') 
ylabel('ffequency [MHz]’) 
figure(2) 
contour(t,fd,z,5) 
axis([-10 10-1 1]) 
xlabel('delay [microsec]') 
ylabelCfrequency [MHz]') 
figure(3) 

plot(t,z((ny+l )/2,:)) 
axis([-30 30 0 12]) 
xlabel('delay [microsec]') 
ylabel('magnitude') 
figure(4) 

plot(fd',z(:,(r+l)/2)>’) 
xlabel('ffequency [MHz]') 
ylabel('magnitude') 
figure(5) 
contour(tfd^,5) 
axis([-1.5 1.5 -0.2 0.2]) 
xlabel('delay [microsec]') 
vlabel('ffequency [MHz]') 
figure(6) 

plot(t^((ny+l)/2,:)) 
axis([-1.5 1.5 0 12]) 
xlabel('delay [microsec]') 
ylabel('magnitude') 
figure(7) 

plot(fd'^(: ? (r+l)/2)') 
xlabel('frequency [MHz]') 
axis([-0.12 0.12 0 6]) 
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% Filename: program6.m 

% Title: Ambiguity function of the step frequency radar by the equation method 
% Date of last revision: 26 Apr 1 996 
% Written by: Paulo A. Soares 

% Comments: This program computes and plots several figures of the ambiguity function of the 
% step frequency radar, using the equation method. 



%- 

% Parameters initialization 



%- 

clear; 

tau=l; 

T=5; 

N=10; 

df=l; 

d=tau/T; 



% clears all variables 
% tau=pulse width 
% T=pulse repetition interval (PRI) 
% N=number of pulses 
% df=frequency step 
% d=duty cycle 



%- 

% Definition of axis and resolution 

% 

nx- 101; % nx=number of points in the time axis 

t=linspace(-N*T,N*T,nx); % time axis definition 

ny=300 1 ; % ny=number of points in the freq. axis 

fd=linspace(-5/T,5/T,ny); % freq.axis definition 

%- 

% Computations 

%- 

p=floor(t/T); 

r=(t/T)-p; 

z=[]; 

form=l:ny 

x=G; 

Md(m); 
for n=l :nx 

PP=P(n); 
rr=r(n); 
fl=f-pp*df; 
f2=f-(2 *pp+rr)*df; 
f3=f-(pp+l)*df; 
f4^-(2*pp+rr+l )*df; 
if rr>=0 & rr<d 
iffl==0 
xl=l; 
else 

xl=(sin(pi*fl *(tau-rr 4: T)))./(pi*fl *(tau-rr*T)); 
end 

ifrem(f2*T,l)=0 

x2=N-abs(pp); 

else 

x2=(sin(pi*f2*(N-abs(pp))*T))./(sin(pi*f2*T)); 

end 

x=[x,(tau-rr*T).*xl *x2]; 
elseif rr<l & rr>=(l-d) 
if 0—0 
xl=l; 
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else 

xl=(sin(pi*D*(tau+(rr-l)*T)))./(pi*f3*(tau+(rr-l)*T)); 

end 

ifrem(f4*T ? l)=0 

x2=N-abs(pp+l); 

else 

x2=(sin(pi*f4*(N-abs(pp+l))*T))./(sin(pi*f4*T)); 

end 

x=[xXtau4^rr-l)*T).*xl*x2]; 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=abs(z); 



%- 

% Plotting of results and axis labeling 

%- 

figure(l) 
mesh(tfd,z),grid 
xlabel('delay [microsec]') 
ylabel('£requency [MHz]') 
figure(2) 
contour(tJfd^,5) 
axis([-6 6 -0.2 0.2]) 
xlabel('delay [microsec]') 
ylabel('frequency [MHz]') 
figure(3) 

plot(t,z((ny+l )/2,:)) 
xlabel('delay [microsec]') 
ylabel(’magnitude') 
figure(4) 

plot(fd',z(:,(nx+l )/2)’) 
xlabel('frequency [MHz]') 
ylabel('magnitude') 
figure(5) 
contour(t/d£>5') 
axis([-0.7 0.7 -0.05 0.05]) 
xlabel('delay [microsec]') 
ylabel(’frequency [MHz]') 
figure(6) 

plot(t^((ny+l)/2,:)) 
axis([-1.5 1.5 0 12]) 
xlabel('delay [microsec]') 
ylabel(’magnitude') 
figure(7) 

plot(fd\z(:,(nx+l)/2)') 
axis([-0.05 0.05 0 0.5]) 
xlabel('frequency [MHz]') 
ylabel('magnitude') 



88 



LIST OF REFERENCES 



1. Merill I. Skolnik, Introduction to Radar Systems, Second Edition, McGraw-Hill, Inc., 
1980. 

2. Donald R. Wehner, High Resolution Radar, Second Edition, Artech House Inc., 
1995. 

3. Byron Edde, Radar. Principles, Technology, Applications, Prentice-Hall, Inc., 1993. 

4. George W. Stimson, Introduction to Airborne Radar, Hughes Aircraft Company, 
1983. 

5. Nicholas C. Currie, Radar Reflectivity Measurement: Techniques & Applications, 
Artech House, Inc., 1989. 

6. James A. Scheer, and James L. Kurtz, Coherent Radar Performance Estimation, 
Artech House Inc., 1993. 

7. Nadav Levanon, Radar Principles, Wiley & Sons, Inc., 1988. 



89 












90 



INITIAL DISTRIBUTION LIST 



1. Defense Technical Information Center 2 

8725 John J. Kingman Rd., STE 0944 

Ft. Belvoir, VA 22060-6218 

2. Dudley Knox Library 2 

Naval Postgraduate School 

411 Dyer Rd. 

Monterey, CA 93943-5101 

3. Chairman, Code PH 1 

Department of Physics 

Naval Postgraduate School 
Monterey, CA 93943-51 17 

4. Professor Gumam S. Gill, Code EC/G1 5 

Department of Electrical and Computer Engineering 

Naval Postgraduate School 
Monterey, CA 93943-5121 

5. Professor Suntharalingam Gnanalingam, Code PH/Gm 1 

Department of Physics 

Naval Postgraduate School 
Monterey, CA 93943-5117 

6. Direc^ao do ServiQo de Forma?ao 1 

Administra?ao Central de Marinha 

Rua do Arsenal 
1100 Lisboa 
Portugal 

7. LT Paulo A. Soares 2 

Rua Ramiro Ferrao 18, 8° Frente 

2800 Almada 
Portugal 



91 



DUDLEY KNOX LIBRARY 
HAVAL POSTGRADUATE SCHOOL 
MONTEREY GA 93843-510 1 



DUDLEY KNOX LIBRARY 




3 2768 00323868 4 



